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Perturbative quark bound states in pNRQCD* 
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In the framework of the QCD effective field theory called potential Non-Relativistic QCD, we explore quark- 
antiquark bound systems that may be dominated by the perturbative interaction and we discuss the extent of 
validity of such a picture. Some phenomenological implications are outlined. 



1. INTRODUCTION 

If we have a chance to understand the QCD 
bound-state dynamics without resorting to the 
lattice or to models jjj, this is likely to be in 
the cases of bound systems whose interaction is 
dominated by the perturbative dynamics. In such 
cases non-perturbative corrections exist but are 
small and can be parameterized in terms of lo- 
cal or non-local condensates. In the following 
I will call such systems 'Coulombic' or 'quasi- 
Coulombic'. These situations are particularly in- 
teresting not only because, as I will show, they 
have phenomenological relevance, but especially 
because they allow us to understand much more 
about QCD. Here, I will exclusively address such 
situations. 

Typically, systems that may be 'Coulombic or 
quasi' are bound states composed only by heavy 
quarks (and gluons). In such cases, at least for 
the lowest states, the characteristic radius (r) (r 
being the qq distance) of the system is quite small 
and this intuitively justifies the idea that they 
are dominated by the perturbative dynamics: the 
system is too small to probe the confinement ef- 
fects, which arise at a scale 1/Aqcd 3> (r). 

Here, I will introduce an effective field the- 
ory called potential Non-Relativistic QCD (pN- 
RQCD) ^,||, that clarifies such picture and al- 
lows us to include systematically perturbative 



and non-perturbative contributions. Since such 
effective theory is constructed to be completely 
equivalent to QCD, all the information that we 
will obtain go to deepen our knowledge of bound 
systems in QCD. For the non-perturbative situa- 
tion see the talk by A. Vairo at this Conference 
land " 
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2. THE SCALES of QUARKONIUM 

As it is apparent from the spectra, heavy 
quarkonia are non-relativistic systems. Thus, 
they may be described in first approximation us- 
ing a Schrodinger equation with a potential in- 
teraction. This amounts to saying that the heavy 
quark bound state is characterized by three en- 
ergy scales, hierarchically ordered by the quark 
velocity v -C 1: the quark mass m (hard scale), 
the momentum mv (soft scale (S)), and the bind- 
ing energy mv 2 (ultrasoft scale (US)). In the 
Coulombic or quasi-Coulombic situation it is v ~ 
a s . 

In perturbative calculations, these scales typ- 
ically get mixed in the Feynman diagrams and 
originate technical complications (the same hap- 
pens in QED, e.g. for positronium). 

In QCD a further conceptual complication 
arises if we take into account the existence of the 
non-perturbative scale Aqcd, at which the non- 
perturbative effects become dominant. For heavy 
quarks only the hard scale m is surely bigger than 
Aqcd and can be treated perturbatively. 

The existence of these different scales makes 
even a purely perturbative definition of the static 
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Figure 1. Scales and EFTs in NR bound state 
systems. 



qq potential not free from complications. Let us 
consider the energy of the static quark sources, 

1 

E s (r) = lim — ln(Wn), (being W a the static 

Wilson loop of size r x T, and the symbol ( ) 
being the average over the gauge fields), which 
is usually considered as a definition of the static 
potential. At three loops E s shows infrared di- 
vergences g^] . These singularities may indeed be 
regulated, upon resummation of a certain class 
of diagrams, which give rise to a sort of dynam- 
ical cut-off provided by the difference between 
the singlet and the octet potential. However, 
such a dynamical scale is of the same order of 
the kinetic energy for quarks of large but finite 
mass and, therefore, should not be included into 
a proper definition of the static potential in the 
sense of the Schrodinger equation. This is similar 
to what happens for the Lamb shift in QED at 
order 1/m 2 . In QCD this effect calls, even in the 
definition of the static potential, for a rigorous 
treatment of the bound-state scales. 

To address the multiscale dynamics of the 
heavy quark bound state, the concept of effec- 
tive field theory turns out to be not only helpful 
but actually necessary. QCD effective field theo- 
ries (EFT) with less and less degrees of freedom, 
can be obtained by systematically integrating out 
the scales above the energy we aim to describe. 
This procedure leads ultimately to a field theory 
derived quantum mechanical description of these 
systems. The corresponding EFT is called pN- 
RQCD §,§(cf.alsof|). Here, all the dynamical 
regimes are organized in a systematic expansion 



and the above mentioned problems are solved. 
In particular we get a procedure to calculate the 
static potential beyond two loops and to obtain 
the energy levels beyond next-to-next-to leading 
order (N 2 LO)0. 

3. NRQCD and pNRQCD 

Since the typical scales of the system are widely 
separated, it is possible to perform an expan- 
sion of one scale in terms of the others and in- 
tegrate out step by step the scales at higher en- 
ergy. Roughly speaking, first one writes the EFT 
as a 1/m expansion and integrates out the scale 
to in QCD: this leads to NRQCD. Later, one 
writes the EFT as an expansion in the inverse 
of the soft scale, the so called multipole expan- 
sion, and integrates out the soft scale in NRQCD: 
this leads to pNRQCD. See Fig. 1. The effective 
theory supplies us with the procedure to make 
this expansion consistent with the ultraviolet be- 
haviour of QCD and consistent with a systematic 
power counting in the small expansion parame- 
ter v. Only if the physical system is such that 
the scales we are integrating out are much bigger 
than Aqcd, namely mv ^ Aqcd, then it is possi- 
ble to take advantage of the further simplification 
of performing the calculations in a perturbative 
expansion in a s and to recover a perturbative de- 
scription of the dynamics. This is the situation 
we consider here. 

The Lagrangian of NRQCD Q can be organized 
in powers of 1/m. In the two-fermions sector it 
is of the type obtained via a Foldy-Wouthuysen 
transformation, but, since we are modifying the 
ultraviolet behaviour of the theory, matching co- 
efficients and new operators have to be added in 
order to mock up the effects of heavy particles 
and high energy modes into the low energy EFT. 
Since the scale of the mass of the heavy quark 
is perturbative, the scale fj, of the matching from 
QCD to NRQCD, mv < \i < m, lies also in the 
perturbative regime. The integration of degrees 
of freedom is done in practice with a matching 



3 The non-relativistic limit described by the Schrodinger 
equation with the static potential is called 'leading order' 
(LO); contributions corresponding to corrections of order 
v n to this limit are called N n LO. LL means 'leading log'. 



procedure i.e. by comparing on shell amplitudes, 
order by order in 1/m and in a s , in QCD and 
in NRQCD. The difference is encoded into the 
matching coefficients that typically depend non- 
analytically on the scale m which has been inte- 
grated out: c ~ Aa s (ln. ^ + B). One works typ- 
ically in dimensional regularization, MS scheme, 
and with quark pole masses. 

After integrating out the soft scale in NRQCD, 
pNRQCD is obtained. The Lagrangian of pN- 
RQCD is organized in powers of 1/m and r 
(multipole expansion). The matching is done 
by comparing appropriate off-shell amplitudes in 
NRQCD and in pNRQCD, order by order in 1/m, 
a s and order by order in the multipole expansion. 
The matching coefficients are non-analytic func- 
tions of r and have typically the following struc- 
ture: V ~ V(r, p, Si, S 2 )(A' In mr+B' In n'r+C). 

4. pNRQCD for mv > Aq C d 

At the scale of the matching \j! (mv ^S> // 3> 
mv 2 , Aqcd) we have still quarks and gluons. We 
denote by R = (xi +X2)/2 the center-of-mass 
of the qq system and by r = xi — X2 the rela- 
tive distance. The effective degrees of freedom 
are: qq states (that can be decomposed into a 
singlet 5(R, r, t) and an octet 0(R, r,t) under 
color transformations) with energy of the order 
of the next relevant scale, 0(Aqcd, mv 2 ), and 
momentum 4 p of order 0(mv), plus ultrasoft glu- 
ons A M (R, t) with energy and momentum of order 
0(AQCD,mv ). Notice that all the gluon fields 
are multipole expanded. The Lagrangian is then 
an expansion in the small quantities p/m, 1/rm 
and in 0(A QC d, mv 2 ) x r. The pNRQCD La- 
grangian is given at the leading order in the mul- 
tipole expansion by||: 
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4 Although, for simplicity, wc describe the matching be- 
tween NRQCD and pNRQCD as integrating out the soft 
scale, the relative momentum p of the quarks is still soft. 
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bound states 


m 


^(25), Xco(lP), - 




T(25), xw(lP), - 




long-range hybrids, ... 


mv 


Vc J/ip, B c ,... 




short-range hybrids, ... 


mv 2 


tt, T(1S),... 


< mv 2 


QED bound states 



Table 1 

An indicative classification of the states which can 
be treated perturbatively up to a scale mv 11 assum- 
ing the criterion m|j™ rt /AQCD ^ 1- 

+.9-^-Tr {O f r • EO + O f Or • E} - jF£ u F^ a . 

All the gauge fields in Eq. (Ill) are evaluated 
in R and t. In particular E = E(R, t) and 
iD O = id O - g[A (R,t),O}. The quantities 
denoted by V are the matching coefficients. We 
call V s and V the singlet and octet static match- 
ing potentials respectively. At the leading or- 
der in the multipole expansion, the singlet sec- 
tor of the Lagrangian gives rise to equations of 
motion of the Schrodinger type. The two last 
lines of (|l|) contain (apart from the Yang-Mills 
Lagrangian) retardation (or non-potential) effects 
that start at the NLO in the multipole expansion. 
At this order the non-potential effects come from 
the singlet-octet and octet-octet interactions me- 
diated by a ultrasoft chromoelectric field. 
Various applications of pNRQCD have been pre- 
sented at previous editions of this Conference P, 
0. 

Recalling that r ~ 1/mv and that the operators 
count like the next relevant scale, 0(mv 2 , Aqcd)> 
to the power of the dimension, it follows that each 
term in the pNRQCD Lagrangian has a definite 
power counting. This feature makes £ p nrqcd a 
suitable tool for bound state calculations: being 
interested in knowing the energy levels up to some 
power v n , we just need to evaluate the contribu- 
tions of this size from the Lagrangian. From the 
power counting e.g., it follows that the interac- 
tion of quarks with ultrasoft gluons is suppressed 
in the Lagrangian by v with respect to the LO ( 
by gv if mv 2 > A QC d). 

In particular, pNRQCD provides us with the 



way of obtaining the matching potentials via the 
matching procedure at any order of the perturba- 
tive expansion ||. In the EFT language the po- 
tential is defined upon the integration of all the 
scales up to the ultras oft scale mv 2 . From the 
matching to NRQCD in the situation Aqcd -C 
mv we can easily obtain the matching potential 
V s at N 3 LL § 



V s (r)=E s (r)\ 



2-loop+N 3 LL 
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where E s is the static energy defined in Sec. 2 (reg- 
ulated by the resummation). We note that V s 
and E s would coincide in QED and that, there- 
fore, this difference here is a genuine QCD fea- 
ture. Such difference is switched on at NLO in 
the multipole expansion. An explicit calculation 
gives f) 
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where /3 n are the coefficients of the beta func- 
tion (a s is in the MS scheme), and a\ and ai 
were given in JO) . We see that the interpretation 
of the potentials as matching coefficients in pN- 
RQCD implies that the Coulomb potential is not 
simply coincident with the static energy E s . The 
Coulomb potential turns out to be sensitive to 
the ultrasoft scale but infrared finite. The same 
happens with the other potentials (like V or the 
potentials that bear corrections in l/m ra ) that can 
equally be calculated via the matching procedure. 
The \i! scale dependence in the potential is can- 
celled for any physical process by the contribution 
of the ultrasoft gluons that are cutoff at the scale 

Then, there are two situations pq[ . If mv 3> 
mv 2 > Aqcd, the system is described up to or- 
der af by a potential which is entirely accessi- 
ble to perturbative QCD. Non-potential effects 
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Figure 2. Graphs contributing to the spectrum up 
to N 3 LL. 

start at order a^ In \j! |l(| . We call Coulombic 
this kind of system. Non-perturbative effects are 
of non-potential type and can be encoded into 
local (a la Voloshin-Leutwyler fl4|| ) or non-local 
condensates: they are suppressed by powers of 
Aqcd/?™^ 2 and Aqcd/wu respectively. If mv 3> 
Aqcd 3> mv 2 , the scale mv can be still integrated 
out perturbatively, giving rise to the Coulomb- 
type potential (ph. Non-perturbative contribu- 
tions to the potential arise when integrating out 
the scale Aqcd H- We call quasi- Coulombic the 
systems where the non-perturbative piece of the 
potential can be considered small with respect to 
the Coulombic one and treated as a perturbation. 
Some levels of tt, the lowest level of bb may be 
considered Coulombic systems H, while the J/ip, 
the rj c and the short-range hybrids may be con- 
sidered quasi-Coulombic. The B c may be in a 
boundary situation, see Table 1. As it is typical 
in an effective theory, only the actual calculation 
may confirm if the initial assumption about the 
physical system was appropriate. 

For all these systems it is relevant to obtain a 
determination of the energy levels as accurate as 
possible in perturbation theory. 

5. CALCULATION of QUARKONIUM 
ENERGIES at 0(ma 5 s In a s ) 

The perturbative energy levels of quarkonium 
are known at 0(ma 4 s ) from @Jl|. pNRQCD 
provides us with a well defined way of calculating 
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and cc ground states have been studied 



the energy levels at higher orders, since the size 
of each term in the pNRQCD Lagrangian is well- 
defined. In order to obtain the leading logs at 
0(ma|?) in the spectrum, V s has to be computed 
at 0(a*ln), vP at 0(a 3 ln), U s (2) at 0(a 2 In) 
and VJ 3) at 0(a s In). The matching to pNRQCD 
at N 3 LL accuracy was calculated in m& and the 
potentials at the requested accuracy were thus ob- 
tained. In Fig.(H) the relevant graphs are shown. 
In particular from the last graph in pNRQCD 
(the singlet turns into an octet double-line by 
emitting an ultrasoft gluon) we extract the de- 
pendence of Vs on /«'• 

The total correction to the energy at 
0(ma 5 s In a s ) is given by the sum of the averaged 
values of the potentials plus the non-potential ef- 
fects, 

5E n>lJ = 8^ ot E n<hj {n') + 6 vs E n>l (ij,') , (3) 



where H = p 2 /2m + V . 

Different possibilities appear depending on the 
relative size of Aqcd with respect to the US scale 
ma 2 . If we consider that Aqcd ~ ma 2 the glu- 
onic correlator in Eq. (g) cannot be computed 
using perturbation theory. We are still able to ob- 
tain all the moc 8 ln(ma s /m) and ma, ln(ma s / //') 
contributions where the p! dependence cancels 
now against the US contributions that have to 
be evaluated non-perturbatively. 

If we consider that ma 2 » Aqcd, Eq. (||) can 
be computed perturbatively. Being ma 2 the next 
relevant scale, the effective role of Eq. (H) will 
be to replace // by ma 2 (up to finite pieces that 
we are neglecting) 
simplifies to 
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The In a s appearing in Eq. (0) come from logs of 
the type In 1/mr. Therefore they can be deduced 
once the dependence on In to is known. The \J 
dependence of Eq. fl) cancels against contribu- 
tions from US energies. At the next-to-leading 
order in the multipole expansion the contribution 
from these scales reads 

5 US E n>l (n') = -*j^T F x (5) 
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plus non-perturbative corrections that can be pa- 
rameterized by local condensates [[L6Ul4|ll| and 
are of order a 2 (^^ )2(Aoc£)2 and higher. For 
the T(IS') and tt non-perturbative contributions 
are expected not to exceed 1004- 150 MeV and 10 
MeV respectively. 

The calculation (p) paves the way to the full 
N 3 LO order analysis of Coulombic systems and is 
relevant at least for tt production and T physics. 
In the first case it is a step forward reaching a 
100 MeV sensitivity on the top quark mass for 
the tt cross section near threshold to be measured 
at future Linear Colliders p3 . In the second case 
they improve our knowledge of the b mass Jl7],p5f . 

In particular, from (m we can estimate the 

N 3 LL correction to the energy level of the T(15) 

and we find 

1730 
SE W1 = — — m b a 4 s (fj,)a s (n')lnl/a s (ii') 

817T 

~ (80 + 100) MeV, (7) 

which appears not to be small. Corrections from 
this N 3 LL terms have been calculated for the tt 
cross section and T(15*) wave function, cf. |l8[ 
and p7|. They also turn out to be sizeable. 



(6) 



Large corrections, however, already show up at 
NLO and N 2 LO and are responsible for the bad 
convergence of the perturbative series in terms of 
the pole mass. This will be discussed in the next 
section. 

6. RENORMALONS, the POLE MASS 
and the PERTURBATIVE EXPAN- 
SION 

The bad convergence of the perturbative ex- 
pansion (for an explicit example see Sec. 7) can 
be, at least in part, attributed to renormalon con- 
tributions. The pole mass, thought an infrared 
safe quantity UM, has long distance contributions 



of order Aqcd|20[|- Also the static potential is 
affected by renormalons (see e.g. p0[). Rephras- 
ing them in the effective field theory language of 
pNRQCD we can say that the Coulomb potential 
suffers from IR renormalons ambiguities with the 
following structure 



V s {r)\ mrcn = C a + C 2 r 2 + 



(8) 



The constant Co ~ Aqcd is known to be 
cancelled by the IR pole mass renormalon 
(2m po ic|iRron = -Co, @). Several mass 
definitions appropriate to explicitly realize this 
renormalon cancellation have been proposed [El| . 
Among the others, the IS* mass[E2| is defined as 
half of the perturbative contribution to the 3 Si 
qq mass. Unlike the pole mass, the IS" mass, con- 
taining, by construction, half of the total static 
energy (2m + v Coul ), is free of ambiguities of 
order Aqcd- Taking e.g. the T(IS'), the IS 
mass is related to the physical T(15) mass by 
E(T(1S)) = 2mis + At- At is the poorly known 
non-perturbative contribution, which is likely, as 
we said, to be less than 100 -j- 150 MeV. In the 
next section we will present an explicit example 
that shows how, using this mass and thus deal- 
ing with quantities that are infrared safe at order 
Aqcd, the pathologies of the perturbative series, 
due to the renormalon ambiguities affecting the 
pole mass, are cured. 

It is possible to show that the second infrared 
renormalon, C 2 — Aqcd' of V s cancels against 
the appropriate pNRQCD UV renormalon in the 
contribution to the potential originating at NLO 



in the multipole expansion. What remains is the 
explicit expression for the operator which absorbs 
the C-2. ~ Aq CD ambiguity (for details see ||] ) . An 
interesting open question is if an explicit renor- 
malon subtraction similar to that one at O(Aqcd) 
can be realized at the subsequent order 0(Aq CD ). 
This may be relevant since this renormalon is re- 
lated to the corrections at N 3 LL discussed in the 
previous section. Indeed, it is still an open prob- 
lem whether the largeness of the N 3 LL corrections 
is an artifact due to our partial knowledge of the 
contributions at this order, or if it is an artifact 
due to the fact that the subsequent renormalon 
cancellation has to be realized at this order or 
finally if it is a true signal of the breakdown of 
the perturbative series. To make more definite 
statements one should know the complete N 3 LO 
or understand the mechanism of cancellation of 
the second renormalon. 

In the next section I will present a concrete 
example of the relevance of the mass renor- 
malon cancellation in order to obtain reliable phe- 
nomcnological predictions. 



7. THE PERTURBATIVE 

TION of the B r MASS 



CALCULA- 



We consider the perturbative calculation up to 
order m<x\ of the energy of the be ground state: 
this will be relevant to a QCD determination of 
the B c mass if this system is Coulombic or at least 
quasi-Coulombic. We also assume the T(15) and 
the J/ijj to be Coulombic or quasi. The question 
if these assumptions correspond to the actual sys- 
tems cannot be settled here. On the other hand 
there is no a priori reason to rule them out. The 
results will tell us how good were our initial as- 
sumptions. 

In order to calculate the B c mass in perturba- 
tion theory up to order af , we only need to con- 
sider the following contributions to the potential: 
the perturbative static potential at two loops, the 
1/m relativistic corrections at one loop, the spin- 
independent 1/m 2 relativistic corrections at tree 
level and the 1/m 3 correction to the kinetic en- 
ergy. We do not consider a s In a s corrections 
since the mechanism responsible for this large 
contributions has not yet been understood. Then, 
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being m rod = m b m c /(m b + m c ), I = 
ln(2CFa s TO rcd //2), J5 (/i) = -m rcd (C F a s (p)) 2 /2 
and p the scale around which we expand a a (r). 

If we use here the pole masses m b = 5 GeV, 
m c — 1.8 GeV and p — 1.6 GeV, then we obtain 
J B(J5 c ) pert ~ 6149 MeV ~ 6800 - 115 - 183 - 353 
MeV, where the second, third and fourth figures 
are the corrections of order a 2 , a^ and a* re- 
spectively. The series turns out to be very badly 
convergent. This reflects also in a strong depen- 
dence on the normalization scale p: at p = 1.2 
GeV we would get E(B c ) peit ~ 5860 MeV, while 
at p = 2.0 GeV we would get E(B c ) peIt ~ 6279 
MeV. The non-convergence of the perturbative 
series (|J) signals the fact that large /?o contri- 
butions (coming from the static potential renor- 
malon) are not summed up and cancelled against 
the pole masses. In order to obtain a well- 
behaved perturbative expansion, we use, now, the 
so-called 15 mass. We consider the perturbative 
contribution (up to order a*) of the 3 5i levels of 
charmonium and bottomonium: 

E(J/i>) peTt = f(m c ); £(T(15)) port = f(m b ), 

which are respectively a function of the c and the 
b pole mass and can be read off from Eq. (||) in the 
equal-mass case, adding to it the spin-spin inter- 
action energy: to(CfO! s ) 4 /3. We invert these re- 
lations in order to obtain the pole masses as a for- 
mal perturbative expansion depending on the 15 
masses. Finally, we insert the expressions m c = 
/-^(J/VOpcrt) andTO h = /- 1 (^(T(15)) pcrt ) in 
Eq. (0). At this point we have the perturbative 
mass of the B c as a function of the J/tp and T(15) 



E(B c ) pc 



Figure 3. 
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perturbative masses. If we identify the perturba- 
tive masses E(J/tp) pert , E(Y(lS)) pci t with the 
physical ones, then the expansion (0) depends 
only on the scale p. The perturbative series turns 
out to be reliable for values of p bigger than 
(1.2 -=- 1.3) GeV and lower than (2.6 -=- 2.8) GeV. 
For instance, E(B c ) pclt = 6278.5 + 35 + 6.5 + 5.5 
MeV at the scale p = 1.6 GeV. Therefore, we 
obtain now a better convergence of the pertur- 
bative expansion and a stable determination of 
the perturbative mass of the B c . This fact seems 
to support the B c being indeed a Coulombic or 
quasi- Coulombic system. By varying p from 1.2 
GeV to 2.0 GeV and A^L =3 from 250 MeV to 350 
MeV and by calculating the maximum variation 
of E{B c ) pcrt in the given range of parameters, we 
get as our final result 



£(B c )pert = 6326±f MeV. 



(10) 



As a consequence of the now obtained good be- 
haviour of the perturbative series in the con- 
sidered range of parameters, the result appears 
stable with respect to variations of p (see Fig. 
|3|) and, therefore, reliable from the perturbative 
point of view. It represents a rather clean predic- 
tion of the lowest mass of the B c . Notice that all 
the existing predictions are based either on po- 
tential models or on lattice evaluation (with still 
large errors). 

Non-perturbative contributions have not been 
taken into account so far. They affect the iden- 
tification of the perturbative masses E(B c ) peTt , 
£(T(15)) P ert) E(J/ip) pei t, with the correspond- 



ing physical ones through Eq.(^). Let us call 
these non-perturbative contributions Ab c , At 
and Aji^p respectively. As discussed before, they 
can be of potential or non-potential nature. In 
the latter case they can be encoded into non- 
local condensates or into local condensates. Non- 
perturbative contributions affect the identifica- 
tion with the physical B c mass roughly by an 
amount ~ — Ajm/2 —At/2 +Ab c - Assuming 



I A 



J/i>\ 



< 300 MeV, | A T | < 100 MeV and A T < 



Ab c < Ajj^p, the identification of our result ( J10J ) 
with the physical B c mass may, in principle, be af- 
fected by uncertainties, due to the unknown non- 
perturbative contributions, as big as ±200 MeV. 
However, the different non-perturbative contribu- 
tions A are correlated, so that we expect, indeed, 
smaller uncertainties. If we assume, for instance, 
At and Aj/^ to have the same sign, which seems 
to be quite reasonable, then the above uncertainty 
reduces to ±100 MeV. This would confirm, in- 
deed, that the effect of the non-perturbative con- 
tributions on the result of Eq. (fiX\ j is not too 
large. 



8. CONCLUSION and OUTLOOK 

We have shown that for heavy quark bound 
systems in the situation Aqcd < mv 2 , the en- 
ergy levels turn out to be calculable in perturba- 
tion theory plus local or non-local condensates. 
pNRQCD provides us with the appropriate tool 
to calculate these energy levels at the desired or- 
der in perturbation theory. The obtained N 3 LL 
contributions turn out to be sizeable. A com- 
plete calculation at N 3 LO is required to settle 
this issue (g6|. For the practical use of the ob- 
tained perturbative series, renormalon cancella- 
tion mechanisms seem to be important. This may 
be relevant also at order a^\na s . In particular, 
I have discussed how a renormalon-free definition 
of the quark masses improves considerably the be- 
haviour of the perturbative series on the concrete 
example of the B c mass calculation. 
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S. Peris (U.A. Barcelona): In your 
Schrodinger equation what definition of the quark 
mass are you supposed to use? 

N. Brambilla: The pole mass. Any other def- 
inition of the mass that does not change the power 
counting can be used, cf. Sees. 6 and 7. 
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